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This article describes experiments performed with an ultracold gas of metastable helium above the 
Bose-Einstein transition. The gas is trapped in an harmonic magnetic potential and collective exci- 
tations are induced. Both the frequency and the damping rate of the lowest monopole-quadrupole 
m = mode of excitation are studied at different temperatures and compared to theoretical pre- 
dictions. Results indicate that the gas goes from a collisionless regime towards a hydrodynamic 
regime as the elastic collision rate increases, when one goes down along the evaporative cooling 
ramp. However we find a discrepancy for the collisional parameters in comparison with predictions 
relying on the value of the scattering length previously estimated. Possible explanations are put 
forward in the final discussion. 



I. INTRODUCTION 

Soon after the experimental realization of Bose- 
Einstein condensation in trapped atomic gases, studies of 
collective excitations in these systems have become a sub- 
ject of intense investigation by several groups At 
very low temperature, when the system is Bose-Einstein 
condensed, it can be described by the hydrodynamic 
equations of a superfluid which give predictions in 
good agreement with experiments performed with alkali 
condensates. At higher temperature, the mean field ef- 
fects become less important and collisions turn out to 
be the dominant factors. For an ultracold gas above the 
critical temperature of condensation, the dynamical be- 
havior of the dilute gas can be described by the Boltz- 
mann equation. In this case two extreme situations can 
occur: - the hydrodynamic regime, where the mean free 
path between the colliding particles is small compared to 
the dimensions of the cold cloud confined in its trap - and 
the collisionless regime where the motion is described by 
the single particle Hamiltonian. In both cases, when ex- 
cited, the system shows well-defined oscillations result- 
ing from the external magnetic confinement, for which 
predictions have been made (see for instance Q or ||] 
and references therein). The frequency and the damp- 
ing rates of the gas oscillations can be calculated for a 
classical gas confined in an harmonic trap, neglecting the 
mean field effects. In reference || results are given for 
the monopole-quadrupole modes of excitation, both for 
the hydrodynamic and for the collisionless regime, with 
an explicit description of the transition between the two 
regimes. Many experimental studies have focused on col- 
lective excitations of such gases. Experiments have stud- 
ied low-lying excitations over large ranges of tempera- 
ture. For all the available data dealing with the thermal 
cloud of trapped alkali gases significantly above the crit- 



ical temperature T c the regime is collisionless, due to the 
relatively low value of the atomic density and of the s- 
wave scattering length ruling the elastic collision rate at 
very low temperatures. Experiments performed at MIT 
with a dense gas of sodium provided oscillations of the 
thermal cloud analogous to first sound and approaching 
the hydrodynamic regime ||. However this regime was 
never completely reached and the complete check of the 
predictions of ||| is still to come. 

When Bose-Einstein condensation of metastable he- 
lium was first achieved in 2001 ||[|, ^ was soon realized 
that the large estimated value of the scattering length (16 
± 8 nm) would favor a situation where the hydrodynamic 
regime could be observed in the thermal cloud above T c . 
In particular our experiment at ENS seemed appropriate 
for such a study because of its relatively large number of 
atoms at the transition (8 x 10 6 ). We thus decided to fo- 
cus on the study of collective excitations in the thermal 
cloud above T c . The present article first briefly recalls 
the experimental set-up and the quantitative results al- 
ready obtained for the relevant parameters: atomic den- 
sity, temperature, elastic collision rate at the end of the 
evaporative cooling ramp. Then we explain the genera- 
tion of excitations in the cigar shaped thermal cloud and 
describe the optical measurements of the oscillations of 
the monopole-quadrupole m = mode generated along 
the weak axis of the cloud. The data for the damping 
and the frequency of the mode are compared with the 
predictions of reference |6| and a discussion follows. In 
particular one finds a discrepancy between the elastic col- 
lision rate estimated from the initial data of |J and those 
actually measured through the present method. 
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II. MAGNETIC TRAPPING AND OPTICAL 
DETECTION 

The production of a Bose-Einstein condensate of 4 He 
atoms in the metastable 2 3 Si state at ENS is first re- 
ported in I! and described in more detail in p0| . The 
metastable atoms are produced from a cryogenic dis- 
charge source at a rate of 10 14 at/s/sr. The atomic beam 
is manipulated by a laser operating on the 2 3 5i — 2 3 P2 
transition of the helium atom at 1.083 /im. After trans- 
verse collimation, deflection and longitudinal decelera- 
tion in a 2 m long Zeeman slower, the beam loads a 
magneto-optical trap using laser beams with a large de- 
tuning compensated by a large intensity, in order to min- 
imise the rate of inelastic Penning collisions induced by 
the laser 

The atomic cloud of typically 10 9 atoms at 1 mK is 
then transfered into a magnetic trap through a multistep 
process Jl0| with an efficiency of about 75%. The mag- 
netic trap that we use for the confinement is sketched in 
figure 0. 




FIG. 1. Magnetostatic trap for the harmonic confinement of the 
cold helium gas. Ql, Q2, C3 are the Ioffe-Pritchard coils, HI and 
H2 are the Helmholtz coils for the bias compression. The detection 
is optical, using the absorption of a probe laser beam imaged on a 
CCD camera. 

Three coils Ql, Q2, C3 of relatively small size realize 
an anisotropic magnetic Ioffe-Pritchard trap. Two addi- 
tional Helmholtz coils HI and H2 reduce the bias field 
in order to increase the radial confinement. A current of 
46.6 A in all five coils produces a 5 G bias field, radial gra- 
dients of 266 G/cm and an axial curvature of 186 G/cm 2 . 
The trap depth is then about 16 mK, the trap frequency 
along the x axis (||) is W||/27r = oj x /2 / k = 115 Hz. The ra- 
dial frequency is Wyjl-K = lo z /2tt = 190 Hz without com- 
pensation of the bias and 988 Hz with a 5 G bias field. 
The cold cloud at thermal equilibrium has a cigar shape, 
elongated along the x axis and with a cylindrical symme- 
try in the radial directions along the y and z axis. When 
required the cloud can be further compressed in the radial 



direction: for this we use another pair of Helmholtz coils 
delivering an additional field along the x axis which re- 
duces the value of the bias field below 5 G down to nearly 
0. The current in the 5 coils can be switched off in 200 /is, 
but eddy currents induced in several pieces of the set-up 
in the vicinity of the cell create relatively strong transient 
field gradients which disappear with a time scale of a few 
ms. One has to take into account this phenomenon which 
can create some difficulties when deducing quantitative 
values from time of flight measurements. 

The optical detection is shown in Fig. [IJ One measures 
the absorption of a laser probe beam passing through the 
atomic cloud with a CCD camera (Hamamatsu C4880) 
with 1.5% efficiency at the wavelength A =1.083 /im. The 
laser source is a DBR diode laser (SDL 60702 HI) oper- 
ated in an extended cavity and amplified by an Yb doped 
amplifier. The probe pulse is long enough (100 /xs) and 
intense enough (/ = 0.21 sat ) to give a reasonable signal 
to noise ratio in spite of the low efficiency of the detector. 
The magnification being of the order of 1, the resolution 
of the images is given by the size of the detector pixel, 
namely 24 /im. Usually the image of the cold cloud in 
the trap is too small to be detected with a good resolu- 
tion. Thus after switching off the trap, we let the cloud 
expand for a few ms and take a picture. We deduce from 
the images both the actual size of the cloud in the trap 
and its temperature. 

The temperature is ramped down in the magnetic trap 
by evaporative cooling. Details about the optimisation 
of the ramp are given in reference ]To| . In 8 s one de- 
creases the temperature down to the /iK range, whereas 
the phase space density increases by 6 orders of mag- 
nitude and one reaches the Bose-Einstein condensation 
around 5 /xK. At the threshold of the transition the di- 
mensions of the cloud have decreased to 90 and 9 /im in 
the axial and radial directions, respectivrly. 

III. RELEVANT PARAMETERS 

A. Number of atoms 

An important parameter for collision studies is the to- 
tal number of atoms. As already discussed in J^] , a proper 
calibration of this number N is difficult to extract with 
a good precision from the absorption images, due to var- 
ious difficulties dealing with the duration of the pulse, 
the Penning collisions that it induces, its pushing of the 
atoms, and stray magnetic fields produced by eddy cur- 
rents. We are currently developing a calibration method 
based on fluorescence measurements that will be pub- 
lished elsewhere. For the time being we prefer to rely on 
the method discussed in Q based on the determination 
of the critical temperature T c for the condensation. T c 
is deduced from measurements of the relative number of 
atoms in the condensate and in the thermal cloud as a 
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function of temperature. The determination of the crit- 
ical temperature does not require knowing the absolute 
number of atoms in each phase. However it can be used 
to evaluate the number of atoms iV c at the transition, 
assuming that one deals with an ideal bosonic gas. If so, 
N c and T c are related by: 



/iK. One easily finds that: 



rr\3 



iV c = \.2{k b T c /hV) 



(1) 



where h is the Planck constant, fcf, the Boltzmann con- 
stant and v the geometrical average of frequencies of 
the trap (482 Hz). Measurements of reference || gave 
T c = 4.7 ± 0.5 ;uK. One thus deduces N c = 8.2 x 10 6 
atoms with an uncertainty of 30%. We neglected the 
possible errors resulting either from mean field interac- 
tions or from quantum correlations (see the discussion of 
note 13 in reference ||). Once calibrated at T c the num- 
ber of atoms N is deduced for any temperature above 



B. Scattering length 

The scattering length a is another important param- 
eter relevant for the present studies, as the scattering 
processes for metastable helium atoms occur only in the 
s-wave channel at very low temperatures. There are sev- 
eral ways of estimating the value of a from the present 
experiment. All of them suppose that the absolute num- 
ber of atoms is known. As the uncertainty on this number 
is large, the value of a is also derived with a large error 
bar. We here chose to rely on the method described in || , 
in which the scattering length is deduced from the size 
of the condensate. Its value is determined from the char- 
acteristic size of the ground state of the trap and from 
the chemical potential /i, which can be calculated from 
the size of the condensate in the Thomas- Fermi limit [[| . 
We then find the value a = 16 ± 8 nm, consistent with 
the experimental measurement of || and theoretical cal- 
culations of JlJ] and jl9| estimating a of the order of 8 
nm (with no error estimate). A recent calculation |20| is 
also compatible with the present experimental value but 
does not give more indication on its precision. 



C. Elastic collision rates 

Finally from the previous values of parameters N c and 
a, one can deduce a value for the elastic collision rate 
T co u in the thermal cloud just at the transition before 
condensation. We use the simple formula: 



r. 



(2) 



where n is the average density and v is the mean relative 
velocity of the atoms at the transition temperature of 4.7 



k h T 



and 



v = 4-\/ — — = 23 cms 



n = 1.4 x 10 13 cm -3 . 



The value of the cross section a is estimated from that 
of the scattering length by: 

a = 8vra 2 = 6.4 x 10~ n cm 2 

which finally results in: 

T coll - 2 x 10V 1 . (3) 

The uncertainty on this value is estimated within a factor 
of 4 which combines those on a, on N c and on T. How- 
ever it can be compared to the low angular frequency 
of the trap (wii = 2tt x 115 Hz), and one finds a ratio 
w\\/T co u = 0.04, much smaller than 1. A collision rate 
higher than the trap frequency also means that the mean 
free path L ~ 1/na of the atoms between collisions is 
smaller than the dimensions of the cloud. Near the crit- 
ical temperature, L ~ 11 /im, which is smaller than the 
size of the cloud along the weak axis a\\ ~ 140 /im. We 
can thus assume that one enters into the hydrodynamic 
regime, an interesting feature for a cold gas above the 
transition, as stated in the introduction. These results 
motivated us to undertake the following studies of the 
collective excitations in the thermal cloud above T r . 



IV. GENERATION AND DETECTION OF 
EXCITATION MODES 

We decided to generate the lowest excitation modes in 
the cold gas. The goal of this study is to observe the 
oscillations of the cold cloud in the axial direction after 
it receives a kick appropriate to generate the monopole- 
quadrupole m = mode. Let us recall that in an 
anisotropic trap the modes I — m — (monopole) and 
I = 2, m = (quadrupole) are coupled and this results 
in two so-called : 'monopole- quadrupole' modes. We study 
the low-frequency monopole-quadrupole mode, for which 
the axial and radial oscillations are in opposite phase. 

The calculations of reference || predict changes for the 
frequency uiq and for the damping Tq of the oscillation 
mode when the collision rate varies from a collisionlcss 
regime to a hydrodynamic regime (see Fig. ||). We thus 
decided to measure the response in time of the elliptic- 
ity of the cloud after the mode is excited. We repeated 
the measurements at various temperatures above T c in 
a large range of collision rates, expecting to reach the 
hydrodynamic regime at the vicinity of the condensation 
transition. 
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The excitation is made by a transient pulse of mag- 
netic field giving a kick to the trapped cloud. Several 
procedures were used and finally we found that the most 
efficient one is the following one: 

- One starts with the gas at rest in the trap at a given 
temperature T above T c in the conditions used for the 
evaporative cooling ramp, namely with a bias field Bq 
of 5 G resulting of an equal current in the 5 coils of the 
trap. 

- One abruptly lowers the bias field Bq along the weak 
axis (x direction) to a fraction of a Gauss by applying an 
appropriate current to the additional bias coils mentioned 
above. This creates a compression in the radial direction, 
the cloud becomes very elongated. The trap is no longer 
harmonic but the atomic density increases dramatically. 

- One simultaneously creates a sinusoidal modulation 
of this compression field with the additional bias coils. 
The amplitude of the short modulation pulse is of order 
0.15 G, its duration of order 15 ms (3 periods). We chose 
a modulation frequency of 180 Hz, intermediate between 
2 vtt and 1.55 which are the expected frequencies for 
the low frequency monopole-quadrupole mode in the two 
extreme cases, the collisionless regime and the hydrody- 
namic regime ||. 

- Just after this pulse one changes the bias field again 
and sets it to about 2 G, a value significantly lower than 
the usual bias field Bq of 5 G. This value is chosen in 
order the obtain the maximum compression of the gas 
compatible with an harmonic trap. On the one hand 
the higher the compression, the larger the collision rate, 
facilitating the entrance into the hydrodynamic regime. 
On the other the trap has to remain harmonic, meaning 
that the condition /jiB/kbT 3> 1 must remain fulfilled, to 
avoid that the intrinsic frequency of the trap is changed 
by the anharmonicity. 

One can show that such a procedure allows to excite 
the low frequency monopole-quadrupole mode m = 0. 

V. EXPERIMENTAL RESULTS, COMPARISON 
WITH THEORY 

We monitor the evolution of the ellipticity of the cloud 
as a function of time after the oscillation is produced. 
Actually the ellipticity oscillates at the same frequency 
as the collective mode with the same damping in case of 
a sufficiently small excitation. A measurement is shown 
in Fig.§ 

From such data one extracts the oscillation frequency 
wq and the damping rate Tq of the excitation mode. 
For instance the results of Fig. g, taken at 3T C , give 
luq/2tt = 191 ± 6Hz and r Q /27r = 25.3 ± 0.8 Hz. This 
experiment is repeated several times at different values 
of the elastic collision rate obtained by varying the final 
frequency of the RF knife of the evaporation ramp. 




t (ms) 

FIG. 2. Oscillations of the ellipticity of the atomic cloud in the 
monopole-quadrupole m = mode. The frequency and damping 
given by the fit are loq/uj^ = 1.66 and Fq/lu^ = 0.22. 

The results are compared with the predictions of the 
classical theory of reference j(| , which are recalled in Fig. 
| 
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FIG. 3. Theoretical predictions of [6] for the low frequency 
monopole-quadrupole mode m = of an elongated cloud of cold 
function of the elastic collision rate. On the upper part 
(Fig. 3a) is plotted the damping of the oscillations, and on the 
lower part (Fig. 3b) the frequency of the oscillations, both in units 
of u)n, the axial (weak) frequency of the trap (Figure taken from 
[6])- 

These results derive both from numerical simulations 
and from analytical calculations, which are in excellent 
agreement. One notes that the ratio ujq/ojh decreases 
from 2 for the collisionless regime to 1.55 for the hydro- 
dynamic regime, whereas Tq/uj^ goes to at both ends 
but reaches a maximum in between the two regimes. The 
results of Fig. 3a and Fig. 3b can be combined to give the 
curve of Fig.||. 
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Fig.|| shows the experimental data we obtained. The 
advantage of such a presentation is that only relative 
values of the measured frequencies are used, which are 
known with good precision. One does not require knowl- 
edge of the collision rates, which is difficult due to the 
uncertainties in the present stage of the experiment. In 
Fig. |I]we plotted the points (called "ENS") correspond- 
ing to the highest collision rates that we were able to 
achieve at a temperature of about 3T C . The comparison 
between experimental data and theoretical calculation al- 
lows us to make an indirect measurement of the elastic 
collision rate in the excited cloud. The highest collision 
rate that we obtained is in the range of 3.5 x 10 3 s _1 . 
One remarks that we were unable to reach the expected 
collision rate of 2 x 10 4 s _1 which is indicated with a cross 
in Fig. ^ and which was estimated independently at the 
transition 0. We also plotted the results obtained by 
the group of W. Ketterle at MIT f§ for cold Na atoms. 
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FIG. 4. Plot of the frequencies and damping rates of the 
lowest monopolc-quadupole mode m=0. The solid line summa- 
rizes the theoretical calculation already mentioned in Figjs| The 
cross corresponds to the estimated value of the collision rate 
T co n = 2 X 10 4 s _1 for metastable helium as derived from [9]. 
The arrow indicates the displacement on the curve when T co u in- 
creases. ENS: Experimental results of this work for metastable 
helium. MIT: Experimental data point of [3] for sodium. 



VI. DISCUSSION 

Usually, when the initial clastic collision rate is high 
enough, the RF evaporation ramp enters in the so-called 
run-away regime, where both phase space density and 
T co u increase exponentially. So we expect to measure 
the highest collision rate at the Bose-Einstein transition. 
In our case, we observed that the phase space density 
always increases (leading to the Bose-Einstein condensa- 
tion) but on the contrary, at the very end of the evap- 
oration ramp, the collision rate starts to decrease. The 
two behaviors are not incompatible since the phase space 



density scales like N/T 3 whereas the collision rate scales 
like N/T. Monitoring the evolution of the elastic colli- 
sion rate along the evaporation ramp can be easily done 
while monitoring the optical density {OD) at the center 
of the cloud. Actually, both quantities are proportional 
in a harmonic trap : OD oc N/T oc Tcoii- The evolution 
of the optical density of the cloud during the evapora- 
tive cooling is plotted in Fig. |[ This plot shows that 
the highest elastic collision rate we can produce is ob- 
tained at about 3Tc- Then, instead of getting deeper 
into the hydrodynamic regime, the elastic collision rate 
turns back towards lower values while approaching the 
transition. 
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FIG. 5. Measurement of the optical density at the center of the 
cloud as a function of the trap depth at the end of the evaporative 
cooling ramp. The optical density being proportional to the elastic 
collision rate in a harmonic trap, this plot shows that one leaves 
the run-away regime below a temperature of about 3T C . 

Actually, while studying the monopole-quadrupole ex- 
citation, the highest elastic collision rate measured in Fig. 
@ (TcM ~ 3.5 x 10 3 s _1 ) was obtained after evaporation 
down to about 3Tc, and after additional compression of 
the magnetic trap. From this value we can infer the cor- 
responding elastic collision rate before additional com- 
pression : 2 x 10 3 s _1 . Then, from Fig. | we can infer 
the collision rate at the critical temperature Tq which is 
about twice as small, leading to 10 3 collisions per second. 
This value is about 20 times smaller than the one derived 
from the measurement of the critical temperature of the 
phase transition. 

Within our error bars, assuming 8 nm for the scattering 
length instead of 16 nm reduces by 4 the discrepancy be- 
tween the two independent measurements of the elastic 
collision rate. Then, when we first estimated the colli- 
sion rate knowing the critical temperature, the number 
of atoms at the transition was nearly twice as large as it 
was when we produced collective excitations, which also 
explains part of the discrepancy. Still it is not enough to 
conclude to the agreement between the two values. We 
are still working on improving the accuracy of both kinds 
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of measurement to understand the discrepancy. In par- 
ticular we started studying in more detail the effects of 
eddy currents on the measurement of the temperature. 
In fact, magnetic fields remaining after the switching off 
of the trap Zeeman shift the optical lines and modify the 
absorption cross section of atoms during the optical de- 
tection. This modification is not homogeneous in space 
and might lead to a systematic error on the effective mea- 
surement of the size of the cloud after free expansion, 
and, consequently, of the critical temperature. 

In order to find out why we leave the run-away regime 
at the end of the evaporative ramp, we studied losses 
and heating rates. For instance we measured heating 
rates of 20 [iK/s at 3T C . A possible explanation is re- 
lated to inelastic Penning Ionization taking place when 
the density becomes high enough. A dominant process is 
the formation of Hc + ions and He atoms in the ground 
state which have very large kinetic energy originated 
from the internal energy of the metastable atoms (about 
20 eV). These hot products might collide with trapped 
metastable atoms and heat them up via elastic collisions 
pl[ . Estimates from numbers given in [|l) lead to the 
right order of magnitude for the heating rate. The more 
hydrodynamic, the more efficient the heating process will 
be J2^] , apparently leading to a limitation in the highest 
elastic collision rate we can produce. 

In conclusion, the route towards the hydrodynamic 
limit has been explored further than in any other re- 
ported experiment on collective excitation modes. Es- 
pecially the measured frequency of the excited mode is 
strongly shifted towards the hydrodynamic regime as pre- 
dicted. The experiment is still under progress, and we 
hope to go deeper into the hydrodynamic regime. How- 
ever, intrinsic limitations due to inelastic collisions have 
been found. 
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